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Abstract. We present an analysis of archival RXTE and BeppoSAX data of the X-ray source 4U2206+54 . For the first time, 
high energy data (> 30keV) are analyzed for this source. The data are well described by comptonization models (CompTT 
and BMC) in which seed photons with temperatures between 1.1 keV and 1.5 keV are comptonized by a hot plasma at 50keV 
thereby producing a hard tail which extends up to, at least, 100 keV. We offer a new method of identification of neutron star 
systems using a temperature - luminosity relation. If a given X-ray source is characterized by a low bolometric luminosity and 
a relatively high color blackbody temperature (> 1 keV) it has necessarily to be a neutron star rather than a black hole. From 
these arguments it is shown that the area of the soft photon source must be small (r « 1 km) and that the the accretion disk, if 
present, must be truncated very far from the compact object. Here we report on the possible existence of a cyclotron line around 
30keV. The presence of a neutron star in the system is strongly favored by the available data. 
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1. Introduction 

High Mass X-ray binaries (HMXRBs) are often subdivided in 
two broad categories depending on the type of the optical com- 
panion: Supergiant systems and Be/X systems. The vast major- 
ity of these High Mass X-ray Binaries harbour X-ray pulsars 
(Bildsten et al. 1997). Around 70% of all known X-ray pulsars 
are in Be/X-ray systems. These are believed to be young neu- 
tron stars (NS) with relatively strong magnetic fields (B ~ 10 12 
G). The X-ray spectra of these objects are usually described by 
powerlaws modified at high energies by a cutoff around 7-10 
keV (Lewin et al. 1995). Some HMXRBs do not display X-ray 
pulsations. From these, three show the typical characteristics of 
black hole candidates, two of them being in the LMC and only 
one (Cyg X- 1 ) in the Galaxy. 

There are, however, other HMXRBs in which pulsations 
have not been detected, in spite of intensive searches, that do 
not display the characteristics of black hole candidates. The 
nature of their compact objects is therefore unknown. One such 
system is 4U2206+54. 

4U2206+54 was first reported by Giacconi et al. (1972) us- 
ing the Uhuru satellite. Steiner et al. (1984) identified the opti- 
cal counterpart with the early type star BD+53°2790. This star 
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showed Ha in emission while the photometric colours where 
consistent with a B0 star and, consequently, the X-ray source 
was proposed as a Be/X-ray binary. It is a persistent X-ray 
source, having been detected by all satellites that have observed 
it. The luminosity is low (Lx ~ 10 35 erg s _I ) and fairly constant 
in the long term but shows flaring with changes of an order of 
magnitude in short timescales. Fits to EXOSAT (Saraswat & 
Apparao 1992) and RXTE (Negueruela & Reig 2001, NR01 
hereinafter) data favored powerlaw models. Using PC A data 
on board RXTE Corbet & Peele (2001, CP01) did not find any 
evidence of pulsations. Instead they found a modulation in the 
ASM X-ray lightcurve with a period of ~9.6 d, which can be 
interpreted as the binary period. 

On the other hand NR01 concluded that the spectral clas- 
sification of the optical companion is far from being well es- 
tablished. The optical spectrum does not belong to any specific 
spectral type and is rather peculiar. These authors conclude that 
BD+53°2790 is not a classical Be star but that the most proba- 
ble scenario is a peculiar active 09p star of moderate luminos- 
ity (III to V) orbited by a compact object which accretes matter 
from the stellar wind. 

The nature of the compact object, in turn, is not clear. The 
presence of a cutoff at ~ 7keV in the X-ray spectrum (NR01, 
CP01) favours a neutron star interpretation. However, the lack 
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Table 1. Journal of the observation for the 1998 November 23 
BeppoSAX observation. Ti nt is the net integration time in ks. 
The net integration time of PDS is shorter by ~5Q% due to 
the high deadtime of the instrument. The net integration time 
of the LECS is ~50% shorter as well due to the fact that it 
observes only during night time. The beginning and the end of 
the observation is in MJD. 



Detector 


c/s 


Start 


Stop 


Tint 


LECS 


0.197+0.005 


51140.710 


51140.853 


12.4 


MECS 


0.404±0.004 


51140.685 


51140.828 


33.4 


PDS 


0.451+0.006 


51140.686 


51140.829 


14.6 



Table 2. RXTE observed 4U2206+54 in 1997 March and 2001 
October. Provided are RXTE observation ID, the beginning 
and the end of the observation in MJD, and the resulting on- 
source time after initial data screening and dead time correc- 
tion for the PCA and the HEXTE in ksec. For details of the 
data extraction, see text. 



ObsID 


Start 


Stop 


PCA 


HEXTE 


20140-01-01-00 


50518.176 


50518.407 


9.1 


5.9 


20140-01-01-01 


50520.251 


50520.365 


4.5 


3.1 


60071-01-01 


52194.010 


52194.164 


17.7 


12.2 


60071-01-02 


52196.268 


52196.601 


29.8 


19.2 


60071-01-03 


52198.834 


52199.017 


20.9 


13.9 


60071-01-04 


52201.339 


52201.433 


19.6 


12.6 


Total 






101.6 


66.9 



of pulsations does not allow to confirm it unambiguously. 
The system displays striking analogies with the microquasar 
RX J1826.2-1450-LS 5039 (optical counterpart 06.5V[(f)], 
low persistent luminosity and no pulsations). In the case of 
4U2206+54 the presence of a black hole would be consistent 
with the above characteristics (especially the lack of pulsa- 
tions). However, as we will show in this paper, there is strong 
evidence for the presence of a NS in the system. 

So far, the spectral analysis has been limited to the 2- 
30keV band. The spectra have been usually fit using powerlaws 
modified at high energies by cutoffs. In order to gain a more 
physical insight into the properties of this mysterious system 
we present in this paper an analysis of partially unpublished 
archival RXTE data and unpublished BeppoSAX data of three 
different epochs. 4U2206+54 was observed in a very broad en- 
ergy range from 0.5 to 200 keV by BeppoSAX and from 2.5 
to lOOkeV by RXTE. For the first time, data from the High 
Energy X-ray Timing Experiment (HEXTE) onboard RXTE 
are analysed for this source. 

2. X-ray observations 

2.1. BeppoSAX data 

We present an observation made in 1998 November 23 
using the Low-Energy Conc entrato r Spectrometer (LECS; 
0.1-10keV; Parmar et al. Jl997l) '). the Medium-Energy 
Concentrator Spectrometer (MECS; 1.8-10keV; Boella et 
al. dl997l) 1 and the Phoswich Detection System (PDS; 15- 



300 keV; Frontera et al. (Il997ln on-board BeppoSAX. All these 
instruments are coaligned and referred to as Narrow Field 
Instruments (NFI). The MECS consists of two (three until May 
9, 1997) grazing incidence telescopes with imaging gas scin- 
tillation proportional counters in their focal planes. The LECS 
uses an identical concentrator system as the MECS, but utilizes 
an ultra-thin entrance window and a driftless configuration to 
extend the low-energy response to 0.1 keV. The non-imaging 
PDS consists of four independent units arranged in pairs each 
having a separate collimator. Each collimator was alternatively 
rocked on-source and 210' off-source every 96 s during the 
observation. The High Pressure Gas Scintillation Proportional 
Counter on-board BeppoSAX (HPGSPC) was not used in this 
observation as the source was not bright enough for this instru- 
ment. Table [2 lists the BeppoSAX observation epoch and net 
exposure times for our source. The LECS instrument is able to 
observe only during night time. Therefore, the on-source time 
for LECS is usually a half that of MECS. 

Good data were selected from intervals when the eleva- 
tion angle above the Earth's limb was > 4° and when the 
instrument configurations were nominal, using the SAXDAS 
2.0.0 data analysis package. In order to produce the spectra 
and lightcurves, LECS and MECS data were extracted cen- 
tered on the position of the source using the standard radii of 
8' and 4', respectively. Background subtraction for the PDS 
was performed in the standard way using data obtained dur- 
ing intervals when the collimators were offset from the source. 
Background subtraction for the imaging instruments (LECS 
and MECS) was performed using blank sky fields provided by 
the BeppoSAX Science Data Center, using the same region of 
the detector as the source. 

The LECS and MECS spectra were rebinned to oversam- 
ple the full width half maximum of the energy resolution by 
a factor 3 and to have additionally a minimum of 20 counts 
per bin to allow use of the x 1 statistic. Data were selected in 
the energy ranges 0.1-10.0keV (LECS), 1.65-10keV (MECS) 
and 1 5-200 keV (PDS) where the instrument responses are 
well determined and sufficient counts are obtained. This gives 
background-subtracted count rates of 0.197, 0.404, and 0.451 
s" 1 for the LECS, MECS and PDS, respectively. 

All uncertainties quoted are on a 90% confidence level for 
one parameter of interest unless otherwise specified. 

2.2. RXTE data 

We used archival RXTE data of 4U2206+54 for this analysis. 
RXTE observed 4U2206+54 in A02 from 1997 March 11 to 
1997 March 12 and four times between 2001 October 12 and 
2001 October 20. The observations are evenly distributed over 
the possible orbit. The 1997 observations were quite short and 
resulted in less then 15 ksec onsource time in the PCA. The 
2001 observation was substantially longer and resulted in al- 
most 90 ksec usable onsource time in the PCA. The HEXTE 
exposure times are substantially shorter due to the rocking of 
the two clusters and the deadtime of the instrument (see below). 
A detailed overview of the observations is given in Table|2] 
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Fig. 1. Left panel.-O.l-lO.OkeV LECS, 1.65-10keV MECS and 15 - 200 keV PDS lightcurves . Data have been rebinned to 60s. 
Right panel: MECS power spectrum. In agreement with previous studies, no significant period is found. The same is true for the 
other instruments. 



RXTE consists of three instruments (for a complete de- 
scription of the satellite, see Bradt et al. 1993): the Proportional 
Counter Array (PC A Jahoda et al 1996), the High Energy X- 
ray Timing Experiment (HEXTE Rothschild et al 1998), and 
the All Sky Monitor (ASM Levine et al. 1996). 

The PCA consists of five co-aligned Xenon proportional 
counter units with a total effective area of ~6000cm 2 and a 
nominal energy range from 2 keV to over 60 keV. However, due 
to response problems above ~20keV and the Xenon-K edge 
around ~30keV, we restricted the use of the PCA to the en- 
ergy range from 3 keV to 20keV (see also Kreykenbohm2002). 
Since the instrument is pointing continuously at the source, 
a background model is used for background subtraction. We 
used the Faint model as is appropriate for a dim source like 
4U2206+54 (see for a description Stark 1997). To account for 
the uncertainties in the PCA response matrix, we used system- 
atic errors. 

On the other hand, we have estimated the galactic diffuse 
background emission at the source position as measured by the 
ROSAT PSPC instrument. Then we converted this count rate 
into fluxes with the PIMMS tool provided by HEASARC. The 
predicted 3-10 keV flux is less than 10~ 14 erg s _1 crrT 2 and 
does not contaminate the source spectrum. 

The HEXTE consists of two clusters of four 
NaI(Tl)/CsI(Na) Phoswich scintillation detectors with a 
total net detector area of 1600 cm 2 . However, early in the 
mission, an electronics failure left detector three in cluster B 
unusable for spectroscopy. These detectors are sensitive from 
~15keV to over 200 keV; however, response matrix, instru- 
ment background and source count rate, limit the energy range 
from 18 to lOOkeV in the case of 4U2206+54. Background 
subtraction is done in the HEXTE by source-background 
rocking of the two clusters every 32 s. 

To improve the statistical significance of the data, we added 
the data of both HEXTE clusters (using a 1:0.75 weighted re- 
sponse matrix to account for the loss of one detector). We fur- 



ther improved the significance of the data by binning several 
channels together. 

2.2.1 . Lightcurves 

In FigH]the BeppoSAX LECS, MECS and PDS lightcurves re- 
binned to intervals of 60 s are shown. Variations of a factor of 
3 can be seen on timescales of ~30min. A search for pulsa- 
tions using Scargle periodgram and Epoch-folding techniques 
(see, e.g., Leahy et al. 1983) failed to result in any significant 
peaks (see Fig^ right panel) in agreement with previous works 
(NR01,CP01). 

The lightcurve clearly presents flaring activity with periods 
of relative quiescence. The erratic flaring seems otherwise to 
be typical of wind accretion, as observed in other HMXRBs 
like Vela X-l (Kreykenbohm et al. 2002, Rreykenbohm et al. 
1999). 

2.2.2. Spectrum 

The overall spectrum of 4U2206+54 was investigated by si- 
multaneously fitting data from the BeppoSAX NFI as well as 
the RXT E PCA and HEXTE data using the xspec package 
(Arnaud i 19961) v 1 1 .2.00). 

As can be seen in Tabled the on-source time is very differ- 
ent for the three SAX instruments. Strictly speaking, the MECS 
is sampling more epochs than LECS, while the source proved 
to be very variable during the observation. However, the error 
bars in LECS data are much larger than those of MECS so that 
this lack of simultaneity has a very small incidence on the fit 
statistics. As a matter of fact, the parameters of the fit are not 
too affected if we reject the LECS data completely. If we reject 
all the MECS data not strictly simultaneous with LECS then we 
obtain larger uncertainties in the fit due to the lack of photons. 
We prefer, therefore, to use all available data. 
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Fig. 2. Left panel: XTE-1997 spectrum fit with a high energy cutoff powerlaw. The parameters are compatible with those obtained 
by NR(01). Right panel: XTE-2001 data with fit with the same model. As can be seen, the spectra is well described by the model 
between 3 and 30 keV. Beyond that point however, the residuals present a large bump. The same is seen in the SAX- 1998 data. 



Table 3. Goodness of fit for additive models during the three 
observations. The;^ 2 is shown as well as the degrees of freedom 
(in parentheses). The name of the models are written in xspec 
terminology. All models include photoelectric absorption. 



Model 


XTE-1997 


SAX- 1998 


XTE-2001 


po 


17(54) 


2.11(115) 


8.45(51) 


highecut(po) 


0.71(52) 


1.32(113) 


1.27(49) 


po + bb 


0.94(52) 


1.24(117) 


0.96(49) 


po + diskbb 


0.98(52) 


1.26(117) 


1.31(49) 


bremss 


3.05(54) 


1.52(115) 


2.59(51) 


raymond 


17.83(54) 


1.87(115) 


10.65(51) 



The photoel ectric absorption cross sections of Morrison 
& McCammon (1983) and the solar abundances of Anders & 
Grevesse ( 1989) were used throughout the analysis. 

Since much of the work has been done so far fitting modi- 
fied powerlaws up to 30 keV the first thing we test is whether 
these models fit properly the spectra if the high energy data is 
taken into account. In Table[3]we quote the^- 2 and degrees of 
freedom for several models. 

In Fig.EJwe present the XTE-1997 and XTE-2001 3-90 keV 
data fit with a high energy cutoff powerlaw. As can be seen, 
the model describes correctly the data between 3 and 30 keV, 
the spectral ranges used in the past. However, the model is un- 
able to describe properly the hard energy tail beyond 30 keV. 
In the SAX- 1998 data, the residuals present exactly the same 
behaviour and, in this case, do not describe properly even the 
low energy bins. 

The analysis yields the following conclusions: 

- The very presence of a hard tail is established for the first 
time, using high energy data from the HEXTE and PDS 
instruments. This hard tail extends at least to 100 keV where 
the signal to noise ratio starts to be poor. 

- Simple models including an absorbed powerlaw and an ab- 
sorbed powerlaw with a high energy cutoff give poor fits 
when the hard tail beyond 30 keV is considered (see Fig0 



Table 4. Parameters for the models of Table [3] with lower xl 
(second and third models) during the three pointings. The H 
column is given in units of 10 22 cirT 2 ). Energies and tempera- 
tures are given in keV. 



Model 


Param 


XTE-1997 


SAX- 1998 


XTE-2001 


highecut(po) 


Nk 


4.5±0.4 


1.1+0.3 


4.6±0.1 




E C h/ 


7.6±0.4 


7.8±0.5 


4.3+0.3 




E/oM 


16.3+1.2 


11±3 


20±2 




r 


1.6+0.1 


1.0±0.2 


1.6+0.1 


po + bb 


N B 


4.4±0.4 


1.0+0.3 


4.8+0.5 




kT 


2.16±0.08 


1.74+0.11 


2.4+0.1 




r 


1.96±0.04 


1.53+0.14 


2.0±0.1 



- Models including a soft blackbody or multicolored disk 
blackbody emission plus a hard powerlaw tail give good 
descriptions of XTE data but not for the SAX- 1998 data 
(FigH right panel, specially at low energies). Note that 
the model with a multicolored blackbody disk results in a 
slightly worse fit (see discussion). 

- Other models with combine thermal models and power- 
laws, result in poor fits. For XTE 2001 data, for instance, 
the Bremss + po model results in a^ 2 of 2.5. 

- Models including only thermal components 
(Bremsstrahlung, Raymond, etc.) are unacceptable. 

Although in some instances the powerlaw + bb model 
could give an acceptable fit, in others do not. In principle, this 
additive model is not physically justified but it is a crude ad- 
ditive approximation of the Comptonization model for which 
a blackbody component and powerlaw component are the 
low energy and high energy asymptotics, respectively (see 
Sunyaev & Titarchuk 1980, hereafter ST80). Thus we use exact 
Comptonization models for model fitting to the data rather than 
any approximation of it. In this work we shall deliberately ap- 
ply models which parameters can be more readily interpreted 
in physical terms and have as few free parameters as possible. 
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Fig. 3. Left panel: SAX-1998 spectrum and folded BMC model. The model is only modified at low energies by photoelectric 
absorption. Right panel: Same data with a powerlaw plus blackbody model. 
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Fig. 4. RXTE-2001 data fit with the BMC model. The model 
parameters are given in Table [5] 



The best fits are obtained with the thermal comptonization 
model (compTT in xspec) (Titarchuk ( 1994), Hua & Titarchuk 
with the addition of a blackbody component at low en- 
ergies , and t he Bulk Motion Comptonization model (Titarchuk 
et al. Jl997l) . Shrader & Titarchuk Jl999l) ). or BMC model. 

This single model reproduces rather well the data both at 
high and at low energies where it must be only modified by 
photoelectric absorption. The best fit parameters for these mod- 
els are given in Table |5] 

Assuming a distance of ~ 3 kpc (NR01), the lu- 
minosity of the source was 7.22xl0 35 erg s~'(XTE-1997, 
'high state'), 5.8xl0 34 erg s _1 (SAX-1998, 'low state') and 
1.78xl0 35 ergs _1 (XTE-2001, 'intermediate state') respec- 
tively. Consistently with earlier observations, no iron line 
around 6.4 keV is detected. 



3. Discussion 

As we have seen, the cutoff powerlaw is no longer a satisfac- 
tory description when data beyond 30keV are taken into ac- 
count. Indeed the spectra show clearly the presence of a hard 



tail which extends to the lOOkeV zone. This hard tail can be 
described by a powerlaw or by comptonization models. 

In contrast, models including only thermal components can 
be discarded (see Table [3j. This fact, along with the presence 
of the hard tail, clearly argue against the accretion onto a WD 
or the interaction between two normal stars as the origin of 
the X-ray emission (a possibility put forward by NR01, ow- 
ing to the unusual optical spectrum). Indeed, colliding wind 
systems present spectra which are rather soft with plasma at 
temperatures of kT ~ 1 keV (Cooke, Fabian & Pringle ( 1987), 
Pollock dl987l) 1. In order to produce measurable X-ray emis- 
sion the very strong stellar winds found in WR stars are re- 
quired. However the optical counterpart of 4U2206+54 is a late 
O star of moderate luminosity (Class III or V). According to 
NR01, the putative companion could be a B star. Such a sys- 
tem can not produce the hard X-ray emission observed in this 
system. The presence of a BH or a NS is, therefore, required. 

First we analyze the parameters deduced from the 'stan- 
dard' thermal comptonization plus blackbody emission model 
(comptt + bb in Table |3. Under this interpretation we first 
note that the temperature of the injected soft photons (kTo) is 
rather high while the luminosity of the source is very small. In 
order to reconcile these two facts, a small emission area must 
be invoked. We therefore compute the radius of the Wien soft 
photon source by equating the bolometric luminosity of the soft 
photon source with that of a black body of area nR\ (In't Zand 
et al. 1999): 



Rv, 



fbol 

0.6«'kp C j— tC&To/I keV)- 2 [km I 



+ y 



(1) 



where y is the comptonization parameter y - 4kT e r 2 /m e c 2 
which gives the relative gain of energy from the inverse 
Compton scattering. f h ° l is the flux of the comptonized com- 
ponent of the model, corrected for absorption, in units of 10" 10 
erg s" 1 cm" 2 . This is obtained by integrating the spectrum. For 
the parameters given in Table|5]we obtain Rw — 1 -98 km (XTE 
1997), R w =* 0.88 km (SAX 1998) and R w 0.96 km (XTE 
2001). 
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Table 5. Parameters for the comptonization models 



Model 



Parameter 



XTE 1997 SAX 1998 XTE2001 



compTT + hb 



Tbb 

kT 
kT e 

t p 

F M 

1 2-10 keV 

F M 

r bol 

CR <fc) 
*?(dof) 



2.04+0.12 
0.76+0.04 
1.50±0.06 
47.79±0.06 
0.50±0.3 

26.8 

67.0 

0.28 
0.66 (49) 



0.30 + 0.15 
0.73+0.05 
1.26±0.07 
47.48±0.06 
0.70±0.3 
4.1 
5.4 
0.43 
1.16(114) 



4.22+0.16 
0.65± 0.04 
1.46+0.06 
46.46±0.06 
0.46±0.3 
12.6 
16.5 
0.34 
0.62 (47) 



BMC 



N H (10 22 cm- 2 ) 0.2±0.1 0.30 ±0.15 



0.2±0.1 



kT col (keV) 

/ 

A N (xlO~ 4 ) 
pW 

1 2-10 keV 

CR <4) 
^(dof) 



1.38±0.09 1.09+ 0.08 1.38±0.09 



1.4±0.23 

»1 
35.3± 0.3 

24.5 

0.40 
0.71 (51) 



0.99+ 0.18 1.26+0.15 



»1 
5.5± 0.2 

4.1 

0.72 
1.14(116) 



»1 
15.7± 0.2 

9.9 
0.47 
0.88 (49) 



<a) unabsorbed flux in units of 10 11 erg s 1 cm 2 
{b) computed as F 20 -60kcv/F 3 -20kcV 
:) unabsorbed flux in units of 10~ 11 erg s -1 cm -2 . The spectral range 
is 2-100 keV for RXTE data and 0.5-100 keV for BeppoSAX data 



On the other hand, we note that the temperature of the 
Compton cloud is very high (kT e ~ 50 keV), of the order of 
those found in BHCs. The source is harder (higher color ra- 
tio) when the luminosity is lower (see table|5}. However, these 
'state transitions', typical for BHCs, have been observed also 
in some pulsars specially in low regimes. 

The best single component model fit is achieved by apply- 
ing the Bulk Motion Comptonization model. This model repro- 
duces rather well the whole spectrum, without the addition of 
other components, both at high and at low energi es. 

The B MC m odel (Titarchuk et al. dl997h . Shrader & 
Titarchuk ill 999h is a general model for comptonization of soft 
photons which uses the Green's (spread) functions for the treat- 
ment of upscattering and which takes the form of a broken pow- 
erlaw. This formalism is valid for any kind of comptonization 
(bulk comptonization in first order (v/c), thermal comptoniza- 
tion in second order (v/c) 2 ) and remains valid up to photon en- 
ergies comparable to the mean plasma energy (m e c 2 ~ 5 1 1 keV 
in the case of bulk motion). This model has been applied suc- 
cessfully to BH systems. The origin of the soft component is 
the innermost disk zone where the gravitational energy of mat- 
ter is released due to viscous dissipation as well as geometric 
compression of the matter. By itself, however, this is not a proof 
for the presence of a BH in the system since it can be applied 
also to accretion onto Neutron Stars. In this case the source of 
soft photons is the surface of the compact object. Either the disk 
or the surface (or both) emit a soft black body like spectrum 
with a characteristic color temperature r co i. The comptonizing 
region (a compton cloud or a boundary layer) must cover ef- 
fectively this zone (i.e. the innermost region of the disk or the 



spot over the surface) in order to be exposed to a high fraction 
of the seed photons. 

One of the three free parameters of BMC model is a power 
law spectral index a related to the Comptonization efficiency. 
When a is smaller the efficiency is higher (for details see 
ST80). A value close to 1 (vg. in SAX-1998 data) indicates 
that the source is undergoing a phase transition from the low- 
hard to high-soft state. These transitions can be caused by the 
redistribution of mass accretion rates between Keplerian (disk) 
and sub-Keplerian components or by an increase of the optical 
depth to for gravitational energy release at the shock (BH case) 
or at the surface (NS case). These spectral transition models 
were studied and discussed in detail by Zel'dovich & Shakura 
(1969), for the NS case, Chakrabarti & Titarchuk (1995) for 
the BH case and Titarchuk & Fiorito (2004), for the BH and 
NS cases. In Figure|5]we show the strong correlation between 
the spectral index a and the color ratio (CR). This plot shows 
that the source is softer when it is brighter (because the upscat- 
tering is less efficient when luminosity increases). That is to 
say, the comptonization parameter decreases when luminosity 
(pressumably the mass accretion rate) increases (Titarchuk & 
Fiorito, 2004). 

The large / value (s>l), which is the ratio of the number 
of photons multiply scattered in the converging inflow to the 
number of photons in the thermal component, indicates that the 
soft spectrum has been fully comptonized and that, as a matter 
of fact, we do not see the source of soft photons at all. The 
Compton cloud completely obscures the BMC region, namely, 
the innermost part of the accretion disk or the surface of the 
NS. 
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10 100 
Channel Energy [keV] 

Fig. 5. Color ratio vs. a spectral index. The numbers represent Fig. 6. a Data and folded model of 4U2206+54 using the 2001- 
the luminosity of the source in units of 10 35 erg s _1 . Note that RXTE observation, b shows the residuals when fitting the data 

with NPEX modified only by photo electric absorption, c same 
as b, but an additional soft black body component is included. 



the source is softer when brighter. 



Again, to reconcile the high color temperature of the soft 
emitting region T co i ~ 1.1 — 1.3 keV with the low intrinsic 
luminosity (Lx ~ 10 35 erg s -1 ) a small emission area must be 
invoked. Indeed, we can make use of the comptonization en- 
hancement factor L/Lo where L is the soft photon source lumi- 
nosity andL is the resulting comptonized luminosity (Titarchuk 
et al. dl997l) . eq. 37). For a ~ 1 this is equivalent to ln(l/xo) 
where xq - 2.8kT C0 \/kT e . Using the BeppoSAX data (a ~ 1) 
from Tableland assuming a Compton cloud of kT e ~ 50 keV 
(deduced from the thermal comptonization model) we can de- 
rive the luminosity of the soft photon source (which we do not 
see) and gives Lq = 1.6 x 10 34 erg s _1 . Such a low luminos- 
ity cannot be produced in a inner accretion disk of this tem- 
perature (r co i). Borozdin et al. (1999) (see Eq. 9 and Fig. 9 
in that reference) show that a disk of temperature ~ 1 keV 
has necessarily a relatively high luminosity of the order of 
10 36 - 10 37 erg s~'fhat contradicts the observable luminosity of 
the soft component. Assuming that this luminosity is radiated 
by a blackbody with an area nR^ and temperature r co i we can 
deduce (following the same argumentation to deduce equation 

1) that R w « 0.6 ^L^{kT coi/l keV)~ 2 [km] where Lq is in units 

of 10 34 erg s _1 . For the values given above this yields 0.7 km. 
This is somewhat smaller than the values deduced previously 
and stresses the need for a small emission area. Clearly, it is 
not possible to have a contribution from a disk in this system 
because the inner region of an accretion disk at these tempera- 
tures would produce a much higher X-ray luminosity. Since the 
radius of the soft photon region is < 2 km, we conclude that the 
only viable candidate is a polar cap on a NS surface. 

The BMC model implies an absorption column which is 
compatible with the expected interstellar value (Nh~ 0.3 X 



10 22 cm 2 , using the relation of Bohlin et al. dl978h and E(B - 
V) st 0.5 deduced from optical data, NR01). This can be ex- 
plained if the resulting comptonized spectrum is not absorbed 
further beyond the comptonizing material close to the X-ray 
source. That is to say, the circumstellar environment (surround- 
ing the optical star) should be tenuous at the site of the compact 
object. This is consistent with the hypothesis that the compact 
object orbits a main sequence 'classical' O star (i.e. with nei- 
ther a strong stellar wind, as in supergiants nor a circumstel- 
lar envelope as in Be type stars) and also with the lack of any 
detectable fluorescence iron line produced in a cold thick sur- 
rounding medium. 

This rarefied medium also contributes to the low luminos- 
ity. Indeed the mass loss ratio between a 09.5V and a 09.51 
star is (De Jager et al. 1988): 



-My 

-Mi 



10 



-2 



(2) 



So, roughly speaking, we have a factor of 10 2 less mate- 
rial to power the X-ray source in a Main Sequence star, like 
BD+53°2790 appears to be, than for a supergiant system like 
Vela X-l and thus also a luminosity lower by the same factor. 

At present, none of the characteristics that can establish un- 
ambiguously the nature of the compact object (pulsations, Type 

1 outbursts, mass function, etc) have been established for this 
system. As discussed before, however, the system lacks an in- 
ner accretion disk, and the emission area is of the order of 1 or 

2 km. This is only consistent with emission from a hot spot on a 
NS surface. Note that this conclusion is based on two different 
models. 
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Fig. 7. Left panel: Unfolded XTE-2001 spectrum of 4U2206+54 showing the possible cyclotron line. We used the BMC model 
to describe the continuum. Left panel: Unfolded SAX-1998 spectrum of 4U2206+54 showing the possible cyclotron line. The 
BMC model is used to describe the continuum. The line is at a slightly larger energy probably due to the coarse binning. 



Other proofs, although circumstantial, also point to the 
NS hypothesis. For instance the negative-positive exponential 
model (NPEX Mihara et al. 1995) gives a comparably satisfac- 
tory fit to the data. This model is commonly used to describe 
the the spectra of pulsars (Mihara et al. 1995, Makishimaet al. 
1999, Kreykenbohm et al. 2002). The NPEX model has been 
successfully fitted to several neutron star spectra like Vela X- 
1 (Kreykenbohm 2002). For 4U2206+54 the resulting fit is 
acceptable (xl= 1-2, see Fig. |6j?) or excellent (xl= 0.6, see 
Fig-IU), using an additional low energy component like a black 
body. The parameters obtained are quite similar to those found 
in pulsar spectra. 

3.1. A cyclotron line at ~ 30 keV? 

One final point deserves attention. In Figs.[5]and0] an absorp- 
tion dip around 30 keV is observed. We were first tempted to 
identify the dip as a cyclotron absorption. The addition of a cy- 
clotron line does, however, not improve the fit. Furthermore, 
it is only one energy bin in both spectra which is below the 
continuum level, so an absorption feature cannot be claimed. 
However, it is intriguing that the same feature can be seen 
in all the spectra, i.e., spectra taken with different instruments 
and different epochs. Furthermore, although the npex model is 
a completely different description of the spectrum, the unex- 
plained "dip" at 30keV is also present (see Fig.[6jl. Therefore, 
we are inclined to think that a very weak absorption feature is 
really present in the data. 

In order to obtain some information about the line we have 
first ignored the energy bins from 25keV to 32keV. Then we 
fit the BMC model. Next we freeze these parameters, add a cy- 
clotron absorption feature (cyclabs in xspec terminology), no- 
tice again energy bins from 25 keV to 32 keV, and refit. Finally 
we thaw all parameters of the model and fit again. The contin- 
uum and line are shown in figure For the XTE-2001 spec- 
trum, the line is centered around Eo = 29 keV, with a width of 
1 keV and an optical depth tq = 0.6. 



Unfortunately the uncertainties are large and the analysis of 
the background does not allow to reject the presence of calibra- 
tion effects so we can go no further on this issue with the avail- 
able data. Further observations of the source at higher spectral 
resolution are required. 

If we interpret this feature as a cyclotron resonant scatter- 
ing feature this would identify the source unambiguously as a 
neutron star with a magnetic field of the order of 3.3 x 10 12 G, 
taking into account the gravitational redshift. As a matter of 
fact we need a strong magnetic field if the material is to be 
funneled to the small pole caps. 

The Alfven radius is inversely proportional to the accre- 
tion luminosity and directly proportional to the magnetic field 
(R M ~ m\ n R 6 2/1 L^ n p 4 ^ 7 where p = BR 3 , Frank et al. 2002). 
4U2206+54 presents a combination of low luminosity and 
thus, presumably, low accretion ratio, and (possibly) a strong 
magnetic field. This radius could then be rather large and help 
to truncate the accretion disk (if any) very far from the compact 
object. 

The lack of the inner accretion disk in the system can ex- 
plain several observational facts: 

- the lack of a fluorescence iron line which would originate 
in the cool parts of an irradiated disk 

- the system presents no radio emission down to 0.04 mJy 
(Ribo et al. priv. comm.). Radio emission would be due to 
radio jets whose formation is strongly hampered without 
the inner disk 

- the lack of outbursts over the years: without an inner disk 
there is no large reservoir of matter which could undergo 
large scale perturbations; instead the source accretes di- 
rectly from the thin steady stellar wind. 

If the system has a strong magnetic field, why do not we 
observe pulsations?. This could be explained by the very hot 
Compton Cloud (kT e ss 50 keV) of this system that can wipe 
out the pulsations (Titarchuk et al. 2002). As a matter of fact, 
we see the radiation fully comptonized while the source of 
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soft photons (presumably the hot spot) is completely obscured. 
Then it is not strange that we miss the pulsations in this sys- 
tem. This leaves open the question, however, why the compton 
cloud si so hot in this system and why we do not observe pul- 
sations in other similar systems. Of course, other possibility is 
that the rotational axis and the magnetic axis might not be off- 
set (or very little offset). Further observations will be necessary 
to asses this issues. 

4. Conclusions 

1. The spectrum of 4U2206+54 can be well described by 
comptonization models. These models show soft radiation 
with temperatures kT ~ [l.l-1.5]keV comptonized by a 
hot plasma at kT ~50keV. This comptonization produces a 
hard tail visible up to, at least, 100 keV. 

2. Application of the BMC Comptonization model to the data 
shows that the soft radiation is completely comptonized. 
Namely the observer does not see the direct blackbody ra- 
diation (presumably) from NS surface. 

3. The temperature of the seed photons, along with the low 
luminosity point to a small emission area (r a 1-2 km 
radius) and contradicts the presence of an accretion disk. 

4. The lack of a fluorescence iron line, the lack of radio emis- 
sion and the lack of outbursts throughout the years also ar- 
gue against the presence of an accretion disk in this system. 

5. Points 3 and 4 put stringent restrictions on the origin of the 
seed photons which must be emitted from a hot spot on a 
NS surface. 

6. The X-ray lightcurve shows no evidence of any coherent 
pulsations at any explored energy range (0.5-60keV). It 
presents the erratic flaring characteristic of accretion from 
the wind of a massive companion. A possible explanation 
for the lack of pulsations is that these are wiped out by the 
very hot Compton cloud. 

7. Finally, an absorption dip around 29keV is present which, 
if interpreted as a CRSF, would identify the system unam- 
biguously as a NS binary. Unfortunately, the significance is 
very low and it is not possible to rule out the possibility of 
calibration effects with the available data. Further observa- 
tions are required to definitely asses this important issue. 
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